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Abstract 



Recently, we introduced an ekpyrotic model based on a single, canonical scalar field that 
OO generates nearly scale invariant curvature fluctuations through a purely "adiabatic mech- 

l/^ anism" in which the background evolution is a dynamical attractor. Despite the starkly 

different physical mechanism for generating fluctuations, the two-point function is identical 
to inflation. In this paper, we further explore this concept, focusing in particular on issues 
of non-gaussianity and quantum corrections. We find that the degeneracy with inflation is 
^ broken at three-point level: for the simplest case of an exponential potential, the three-point 

^ amplitude is strongly scale dependent, resulting in a breakdown of perturbation theory on 

small scales. However, we show that the perturbative breakdown can be circumvented — 



and all issues raised in Linde et al. (arXiv: 0912. 0944) can be addressed — by altering 



the potential such that power is suppressed on small scales. The resulting range of nearly 
scale invariant, gaussian modes can be as much as twelve e-folds, enough to span the scales 
probed by microwave background and large scale structure observations. On smaller scales, 
the spectrum is not scale invariant but is observationally acceptable. 
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1 Introduction 



There are two known cosmological phases that transform inhomogeneous and anisotropic initial 
conditions into a smooth and flat universe, in agreement with observations. The first is inflation [1], 
a period of accelerated expansion occurring shortly after the big bang, which requires a component 
with equation of state w < —1/3. Alternatively, flatness and homogeneity can be achieved during 
an ekpyrotic phase [21 El IH El El El , a period of ultra-slow contraction before the big bang, driven 
by a stiff fluid with w > 1; the ekpyrotic phase also suppresses chaotic mixmaster behavior [HI |9l 
[TOl [TT] . See [T2l [T3] for reviews. In both cases, phases with nearly constant w can be achieved 
with a single canonical scalar field with appropriately chosen potential V{(f)). Whereas inflation 
requires V to be flat and positive, an ekpyrotic phase occurs for an exponentially steep, negative 
potential. A fiducial example is a negative exponential potential, V^(0) = —Voe~'^'^^^^^\ with 1, 
corresponding to w = c^/2 1. 

With nearly constant w and a single scalar field, inflation also generates a nearly scale invariant and 
gaussian spectrum for (, the curvature perturbation on comoving hypersurfaces [HI [151 [16] . The 
same is not true for ekpyrosis. Although in Newtonian gauge the scalar field and gravitational po- 
tential fluctuations are scale invariant, these project out of ( [4J. Barring some higher-dimensional 
or stringy effects near the bounce that mixes gravitational potential and curvature fluctuations 
[6] , the resulting spectrum for ( has a strong blue tilt [IHl [T7] . A scale invariant spectrum can be 
obtained with two ekpyrotic scalar fields, through an "entropic mechanism" [THl UHl 120] that first 
produces entropy fluctuations and then converts them to ( [ISl Ell 1221 [23 [211 [23 [2S] , as in the 
New Ekpyrotic scenario [2n[22|l23]. 

Another desirable property of an inflation phase with a single scalar field is that it is a dynamical 
attractor. On super-horizon scales, ( measures differences in the expansion history of distant 
Hubble patches [16]. Because ( ~ 6a /a — )■ constant at long wavelengths in single-field inflation, 
the perturbation can be absorbed locally by a spatial diffeomorphism ^2ZJ- In other words, the 
background solution a{t) is an attractor. 

Achieving scale invariance and dynamical attractor behavior in alternative scenarios has proven to 
be challenging. A contracting, dust-dominated universe yields an equation for ( that is identical to 
inflation [281 [291 [30] ; but ( grows outside the horizon, indicating an unstable background. Similarly 
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for contracting mechanisms that rely on a time-dependent sound speed The contracting 

phase in the original ekpyrotic scenario [H El S] is an attractor |H1 [ID] , but the resulting spectrum 
is strongly blue rather than scale invariant pTl ITO] . In the entropic mechanism relying on two 
ekpyrotic scalar fields |18|, [191 120] , the spectrum is unambiguously scale invariant, but the entropy 
direction is tachyonically unstable ^3] . 

Recently we have proposed a counterexample, the adiabatic ekpyrotic mechanism |32], in which 
a single canonical scalar field drives a contracting background that is a dynamical attractor and 
generates a scale invariant spectrum for (. The mechanism relies on relaxing the usual assumption 
that the equation of state parameter e = —H/H^ = 3(1 + w)/2 is nearly constant, and obtains 
for fairly simple forms of the potential, such as a lifted exponential: 

V{<P) = Vo{l- e-^<^/^^") , (1) 

with Vq > and c ^ 1. The regime of interest is the transition when e rises rapidly from 
e <^ 1, where the constant term dominates, to e ~ ^ 1, where the negative exponential term 
dominates. During this transition, the scale factor is nearly constant, while the equation of state 
parameter varies rapidly as e ~ 1/r^, where r is conformal time. The quantity z = a(r)A/2e(r), 
which determines the evolution of (, therefore scales as z ~ — exactly as in infiation, 

where e ~ const, and a(r) ^ l/(— r)! The two-point function is, therefore, identical to infiation. 

In fact, a recent study has shown that the only single-field cosmologies with unit sound speed 
that generate a scale invariant spectrum for ( and are dynamical attractors consist of [33l l3l] : 
i) infiation, with a(r) ~ Vkl ^"^^ ^ ~ constant; ii) the adiabatic ekpyrotic mechanism [32] 
mentioned above, with e ~ 1/r^ on a slowly contracting background; and in) the adiabatic 
ekpyrotic mechanism on a background that first slowly expands, then slowly contracts. Here we 
focus on the contracting version of the adiabatic ekpyrotic phase; its expanding counterpart will 
be studied in detail elsewhere 05]. See [36] for related work. (Another counterexample proposed 
recently relies on a rapidly- varying, superluminal sound speed Cs(r) [33 EEl ES]- See [211 HO] for 
earlier related work. Even though the background is non-infiationary, ( is amplified because the 
sound horizon is shrinking.) 

The aim of this paper is to further explore the phenomenological implications of the adiabatic 
ekpyrotic mechanism, focusing in particular on non-gaussianities and strong coupling. We show 
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that the degeneracy with inflation is broken by the three-point function. Unhke the highly gaussian 
spectrum of single-field slow-roll inflation, the rapidly- varying equation of state in our case results 
in a three-point amplitude that is strongly scale dependent and peaks on small scales. For the 
simplest potential Q, in particular, non-gaussianities are of C(l) on the largest scales and grow 
as k"^, resulting in a breakdown of perturbation theory on small scales. Moreover, loop corrections 
dominate over the classical computation on small scales, indicating strong coupling. 

However, these pathologies all result from maintaining the transition phase with large c longer 
than necessary. Strong coupling can be circumvented by modifying the potential such that the 
exponent decreases smoothly from c to 6 ^ c after the transition phase has already generated 
an acceptable range of scale invariant fluctuations. We find that as a result the power spectrum 
acquires a strong red tilt on small scales and the two-point amplitude is exponentially small. 
Suppressing the small-scale amplitude in this way enables perturbation theory to be valid on all 
scales, both classically and quantum mechanically. By the same token, this modification also 
addresses all criticisms raised by [41j. Satisfying all requirements comes at a cost, though: the 
range of nearly scale invariant and gaussian modes is now limited, spanning at most a factor of 
10^ in k space, or a dozen e-folds, which is sufficient to account for microwave background and 
large scale structure observations. 

The paper is organized as follows. We begin in Sec. [2] by reviewing the background dynamics 
for the adiabatic ekpyrotic mechanism. In Sec. |3j we calculate the two-point function, confirming 
that the power spectrum is scale invariant, and briefly discuss the scalar spectral index and 
the tensor spectrum. In Sec. |4| we establish that the background is a dynamical attractor by 
showing that various physical observables become smaller in time and approach the background 
solution. Section [5] focuses on non-linearities and higher-order correlation functions. We explicitly 
compute the three-point amplitude for the fiducial potential ([T]) and find that it is strongly scale 



dependent (Sec. 5.1). Although most contributions can be well approximated by the horizon- 



crossing approximation, surprisingly this method fails for the dominant contribution, which instead 



peaks at late times (Sec. 5.2). This growth in non-linearities results in a breakdown of classical 



perturbation theory (Sec. 5.3) and quantum strong coupling (Sec. 5.4). In Sec. [6] we show how 
these problems can be avoided by modifying the potential as mentioned above, derive various 
parameter constraints to ensure that non-gaussianities and quantum effects are under control, 
and offer a few working examples. We briefly review our main results and discuss prospects for 
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future directions in Sec. [7l 



2 Background dynamics 

The adiabatic mechanism is most simply reahzed with the hfted exponential potential ([T]), where 
Vo > and c ^ 1. This potential is approximately constant and positive at large positive 0, and 
tends to a negative exponential form for large negative (f). The example is not designed to represent 
a complete cosmological model; rather, we focus only on a particular range of (j) to illustrate the 
basic idea behind the adiabatic mechanism. The regime of interest is the transition when the 
equation of state rises rapidly from e ^ 1, where the constant term dominates, to e ~ ^ 1, 
where the negative exponential term dominates. We refer to this as the transition phase. Similar 
behavior over this narrow range of A0 ^ Mpi can be obtained for a wide range of potential 
functions V{(j))] we will see that this freedom enables ways of avoiding problems encountered if 
this first example is considered for all (f). 

The form for l^(</)) reminds one of inflationary examples, but our mechanism is emphatically not 
inflationary in nature. This can be seen in different ways: i) the universe is slowly contracting 
just prior to the onset of the adiabatic mechanism; ii) scale invariant perturbations are generated 
within one Hubble time, hence the universe is essentially static in the process; and Hi) because 
e is changing rapidly during the transition, our background evolution violates the usual slow-roll 
condition ^ 1 of inflation. 

Because the scalar field is falling off a steep potential, the evolution is insensitive to the slowly- 
contracting background and is, therefore, driven by the potential: 



In particular, the evolution of is oblivious to the constant term Vq. As a consistency check, we 
will see shortly that the transition phase occurs in less than one Hubble time. Assuming negligible 
initial kinetic energy, the solution is of the standard ekpyrotic form O |3], Hj 



c 



(2) 



Mpi 




(3) 
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(Here — oo < t < 0, with t = corresponding to — )■ — oo.) 



The cosmological background, meanwhile, can be inferred from the H equation 

1 :o 2 

with solution 



HuUt) = ffo + |j ■ (5) 

At sufficiently early times, H is nearly constant, with the constant Hq fixed by the Friedmann 
equation: SH'^Mp^ = 0^/2 + V{(f)) ~ Vq. During the transition phase, the scalar kinetic energy 
nearly cancels the negative exponential term in the potential, leaving the constant term Vq as the 
dominant contribution to H: 

As the universe contracts and |t| decreases, eventually the constant term no longer dominates, and 
the solution approaches a standard ekpyrotic phase with H{t) ~ 2/c'^t. The end of the transition 
phase — and the onset of the ekpyrotic scaling phase — occurs when the time-dependent and 
constant terms become comparable, that is, at 

^end— trail = ^beg— ek n i tt i " \ / 

C l-nol 

The transition phase is also finite in the past. The above derivation neglects gravity, which is a 
poor approximation for sufficiently large positive where the potential is fiat and Hubble damping 
is important. Specifically, the approximation H(j) <^ cVqC''^'^^'^^^ /Mpi implicit in ^ is consistent 
as long as 

t ^ ^beg— tran = TTFT ' (^) 

The transition phase, defined by tbeg-tran < t < tend-tran, therefore lasts less than a Hubble time, 
as claimed earlier. The scale factor, 

atran(t) ^ 1 + Hot + ^ log (Hot) , (9) 

is approximately constant throughout, and the universe is nearly static. (In integrating ^ to solve 
for a(t), we have chosen the integration constant such that the log term vanishes at tbeg-tran-) 
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Given Q, ^ and ([T]), the equation of state parameter can be expressed as 

H 2 1 



(10) 



tran/ 



Deep in the transition phase, |t| ^ l^end-tranl; the equation of state is rapidly varying, e ~ 1/t^ — 
the key to generating a scale invariant spectrum of curvature perturbations. For |t| ^ |4nd-tran|, 
meanwhile, the equation of state tends to a large constant, e — )■ ^ 1, consistent with an 
ekpyrotic scaling phase [21 H] : 

2 

ascaling(t) ~ i-tf^"' ; i^scaling(t) ~ ^ • (H) 

Over the course of the transition phase, therefore, e grows by a large factor: from (9(l/c^) to 

The adiabatic mechanism relies on an exponentially growing e{t) and nearly constant a{t), the 
exact opposite of the exponentially growing a{t) and nearly constant e characteristic of inflationary 
cosmology. In particular, the rate of change of e is never small: 

1 dine 2 ^cnd-tran^ 



H dt (t + tcnd- 



,2 
tran ^ 



ranges from 0{1) to 0{(?) ^ 1 during the transition, so the usual slow-roll condition r] <^ 1 is 
violated throughout. 



3 Power Spectrum 



The generation of perturbations is conveniently described by C, [IHIIS], the curvature perturbation 
in comoving gauge, 50 = 0, hij = a'^{t){l + 2Q6ij, which completely fixes the gauge. The quadratic 
action governing ( is 



C - (VC)^ 



(13) 



where primes denote derivatives with respect to conformal time r, and z is defined as usual by 



z = a{T)^/2^ 
6 



(14) 



It is convenient to work in terms of the canonically-normalized variable, v = zC^^ whose mode 
functions Vk satisfy 

v'k + (e --]vk = Q. (15) 



In the "transition phase", the scale factor is nearly constant, a(r) ~ 1, — this is the slowly 
contracting background typical of ekpyrotic cosmology — hence conformal time and cosmological 



time are approximately the same: t ^ t. It follows that z ^ ^J2e{t), with e(t) given by (10), and 
therefore 

V, + (e - ——^ ^ v, = 0. (16) 

\ v'' ~r ''end— tran J / 

This equation is valid throughout the transition phase and deep in the ekpyrotic scaling phase. 
The approximation a(t) ^ 1 eventually breaks down as t — )■ 0~, since a(t) ~ (— t)^/*^^, e = 
and hence z/z = in the ekpyrotic scaling phase. Comparison with z/z = 2/tl^^_^^^^ 



implied by the late-time limit of (16) shows that our mode function equation breaks down at 



t ~ tend-tran/c- But this is Well after all modes of interest have exited the Hubble horizon. 



Therefore, (16) accurately describes the generation and freeze-out of scale invariant modes during 



the transition phase, as well as their Hubble exit during the ekpyrotic scaling phase. 



For \t\ ^ |tond-tran|, ( pl6| ) reduces to the same mode function equation as in inflation, where e is 
nearly constant and a{t) grows exponentially! The two-point function thus generated is therefore 
identical to inflation and, in particular, is scale invariant. Indeed, assuming the usual adiabatic 
vacuum, the solution for the mode functions is 

Vk{t) = ( 1 - ^ . (17) 

—trail/ / 

Translating back to the curvature perturbation, (j^ = Vk/z, we obtain 

^'^'Cfc = [1 + ^^(^ + ^end-tran)] e'^'^* . (18) 

2v 2iVfpi 

The corresponding power spectrum as t — )■ 0, defined by (CfcCfc') = (27r)^5^(/c + k')27i'^k^^P(^{k), is 

Pdk) = (1 + ^'4d-tran) • (19) 

The spectrum is therefore nearly scale invariant for /i;|tend-tran| ^ 1, corresponding to modes that 
freeze out during the transition phase. For A;|tend-tran| ^ 1, corresponding to modes that freeze out 
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during the ekpyrotic scaling phase, the spectrum is far from scale invariant, ~ k"^, consistent 
with the strong blue tilt for ( of the original ekpyrotic generation mechanism [H [SJ [TTl [TU] . (The 
blue tilt can be modified by choosing a different V{(f)), as described below.) The range of scale 
invariant modes is determined by the duration of the transition phase: 



^beg-tran _ ^2g^ 



As we will see shortly, requiring that the scale invariant range overlaps with the largest observable 
scales today and that the magnitude of ( matches observations force c to be exponentially large. 

Our analytical treatment is borne out by numerical analysis. Using z = c 

cover the transition and scaling phases, we integrate (15) with c = 140 and \Hq\ = 10~^ over the 



interval —5 x l-f^oT^ < t < —10"^ x \Hq\~-^, over the range of modes 10~^ x \Hq\ < k < 10^ x \Ho\. 
Figure [T] shows the resulting spectrum. The shortest-wavelength modes are barely outside the 
Hubble radius by the end of the integration, which occurs deep in the ekpyrotic scaling phase. 
Modes with k < \Hq\ = 10~^ begin outside the ^-horizon at the initial time and hence are not scale 
invariant. The numerical results show good agreement with the expected range 10~^ ^ ^ ^ 10 of 
scale invariant modes. 



The mode function solution (18) clearly tends to a constant at late times, ( — )■ constant as t — )■ 0, 
again as in inflation. Since ( represents a perturbation of the scale factor in this limit [16], ( ~ 
6a/ a, this implies that the "transition phase" evolution is a dynamical attractor. This statement 
will be made more precise in Sec. |4j but, for the moment, let us contrast this with a contracting, 
dust-dominated universe, corresponding to a(r) ~ (— t)^ and e = 3/2. This background has often 
been hailed as the dual to the inflationary mechanism [HI EHl [29] since z"/z = 2/r^ in this case, 
exactly as in inflation. But because z ~ (— t)^, the curvature perturbation grows at late times, 
C ~ l/(— r)^, indicating that the background is unstable. By contrast, the spectrum generated 
by a slowly-contracting universe with rapidly-varying equation of state, as proposed here, has 
identical two-point function and long-wavelength evolution for ( as inflationary cosmology. 



to 
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Figure 1: Numerical computation of the perturbation amplitude vs. k generated by the 

adiabatic mechanism. The behavior of modes with larger and smaller k depends on the larger 
scenario in which the mechanism is embedded and beyond the consideration of this paper. The 
range of scale invariant modes is in good agreement with the analytical treatment. 

3.1 Observational Constraints 

Our power spectrum must meet two observational criteria. Firstly, the amplitude of the power 
spectrum over the scale invariant range (fc|tend-tran| ^ 1) must match the observed value 



167r2M2j 



10 



-10 



(21) 



which, given Hq, fixes c. Secondly, the scale invariant range must overlap with the scales probed 
by cosmic microwave background and large-scale structure observations, which requires that the 
comoving scale l/|-ffo| encompass the entire observable universe. Assuming that the magnitude of 
H at the end of the ekpyrotic phase, |ifek-cnd|, is comparable to that at the onset of the expanding, 
radiation-dominated phase, then we must demand 



|-^ek— end| ^ '^ek— end | -^ck— end | 

^ 

l^ol ~ '-^today -^today 
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\H, 



ek— end 



today 



(22) 



where in the last step we have assumed a radiation-dominated evolution until the present epoch, 
for simplicity. In other words, 



-f^ck-cndl-f^today ~ 10 \/ \ if gk-end | ^Pl 5 (23) 

which constrains Hq in terms of the reheating scale. For Grand-Unified (GUT) reheating scale, 
|ifek-end| ~ 10^^ GcV, the abovc condition is satisfied for \Hq\ ~ 10"'^ meV, corresponding to 



Vq ~ (10 GeV)^. The constraint (21) on the power spectrum amplitude then fixes c = 10^®. For 
electroweak (EW) reheating, |ifek-end| ~ meV, we similarly get Vq ~ (10~^ MeV)'* and c = 10*°. 

Thus our mechanism requires exponentially large values of c. To recap, this is because the universe 
is nearly static during the generation of perturbations, hence the Hubble parameter must be small 
relative to Mpi for perturbations to overlap with the scales probed by observations. This in turn 
requires a very steep exponential potential in order for the amplitude to match observations. 
Although we have focused on pure exponential potentials, for simplicity, the exponentially large 
values of c required could be achieved, for instance, in the Conlon-Quevedo potential [12], V{(j)) ~ 
exp(— 70*/^), for large 0. 



3.2 Other Observables 

To conclude this Section, we briefly comment on two other observables, namely the scalar spectral 
index and the tensor spectrum. Since the values of c of interest are exponentially large, the 
departures from scale invariance are unobservable for the potential considered thus far. However, 
in this overly simple example, the ekpyrotic phase never ends and the universe does not bounce. 
In a complete model, the exponent c is replaced by c(0) which is exponentially large during most 
of the ekpyrotic phase but is made to fall below one at some given in order to end the ekpyrotic 
phase. The variation in c(0) results a small red tilt, as favored by observations 

V{(j)) = Vo(l- e-/<i<^^W/^^pi) . (24) 

As shown in (22], the resulting tilt is 

n,-l = -4Mpi^. (25) 
10 



Since is decreasing in our solution, the spectral tilt will be slightly red if c^^ > 0. For instance, 
if c(0) changes smoothly by 0{c) during the transition, then 

CA0 log(tcnd-tranAbeg-tran) log(c2/2) ' 

which gives Ug ~ 0.98 for 10^° > c > 10^*, ranging from EW to GUT-scale reheating. Allowing 
for various 0{1) factors that were dropped in this estimate, the generic answer is 1 — ~ few %, 
in good agreement with observations of the cosmic microwave background ^3] . 



Gravitational waves, meanwhile, are not appreciably excited because the background is slowly 
evolving. Tensor perturbations maintain their adiabatic vacuum normalization, hk ~ l/\/k, 
resulting in a strong blue tilt for their spectrum 

Ph{k) ^ ^ k\ (27) 

corresponding to nx = 2. As in earlier renditions of ekpyrotic cosmology |2l HI], the primordial 
tensor amplitude is therefore exponentially suppressed on the largest scales. The dominant grav- 
itational wave background at long wavelengths is the secondary gravitational waves induced by 
the energy density fluctuations, roughly 10~^ times smaller than the primordial fluctuations |45] . 
Detection of primordial gravitational waves, for instance through cosmic microwave background 
B-mode polarization, would unequivocally rule out this mechanism. 



4 Stability 



The fact that ( — )■ constant at long wavelengths suggests that our background is a dynamical 
attractor [27] . To make this statement precise, we will show below that all physical observables 
become smaller in time and approach the background solution. We focus on the transition phase, 
since the ekpyrotic scaling phase has already been shown to be an attractor [H [9l |TOl [TT] . 



Following [lOj, we find that synchronous gauge, in which qqq = —1 and got = 0, is a well-suited 
coordinate system to study stability. The scalar perturbations in this gauge are encoded in the 
spatial components of the metric. 



9ij 



1 + 2C - 2if 



6ij + 2didj 



dt' 



V2 am 



1 



dt' 



H 



(28) 
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as well as in scalar field fluctuations: 

6(f) = 



(29) 



We first show that all metric components tend to their unperturbed value, up to rescaling of 
coordinates, while scalar field perturbations tend to zero. First note that the growing mode 



solution to (15) has the following long-wavelength expansion 

Ck = Co (i + Ikh^ + Oikh')^ , (30) 

where the time-dependent amplitude Co is determined by initial conditions. For the matter per- 
turbation 6(j), this implies 

S<P^^,k'Cot, (31) 



which becomes increasingly small in time. For the metric, the coefficient of the 6ij term in (28) 
gives 

1 + 2C - 2if j dt'j^ ^ 1 + 2Co + . . . , (32) 

where the ellipses indicate terms that become negligible in time. Thus this term goes to a constant, 
which can be absorbed in a (time-independent) spatial gauge transformation. 



Furthermore, since e ~ 6M|j/c^Vot^ during the transition phase, the didj term in (28) becomes 

2A;,A;,Co 2 



2didj J dt' 



log (Hot) + ... , (33) 



where a suitable spatial diffeomorphism has been assumed to normalize the log. The growth of 
this log term looks at first sight dangerous, but note that its time-dependence exactly matches 
the log term in a{t) — see ([9]). This contribution therefore renormalizes the departure from de 
Sitter space of the background solution, and, as such, does not signal an instability. 

It is also instructive to study the stability of curvature invariants. Starting with the Ricci curvature 
of the 3-metric, only the 6ij term contributes to this quantity since the d^dj term is pure gauge: 

r^-(l + 2C-2H f dt'i-\ -> ^Co . (34) 
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The 3-curvature thus goes to a constant, which is acceptable. (The amphfication to a constant 
/^(3) jg precisely how we generate scale invariant perturbations.) Similarly, the perturbation in the 
extrinsic curvature tensor, Kij = gij/2, is of order 

6K 



K HI 



+ . . . , (35) 



where we have used K = 3H{t) = 3{Hq + 2/c^t). Thus, the perturbation in K thus also tends to 
a constant at late times. 

These results, together with earher analyses of standard ekpyrotic scenarios [8|, HI [TOl [11] , establish 
that our cosmological background, consisting of a transition phase followed by an ekpyrotic scaling 
phase, is an attractor solution. The breadth of its basin of attraction is a question that requires 
numerical analysis and will be discussed elsewhere |35j. The fact that small scale modes are 
highly non-gaussian and that non-linearities grow after freeze-out, as discussed in the next Section, 
suggests that the basin of attraction is limited to small perturbations. Moreover, we will see in 
Sec. |6] that our mechanism can produce at most a dozen e-folds of nearly scale invariant and 
gaussian modes. This range, while sufficiently broad to account for large scale observations, does 
not leave much room to wash out arbitrary initial conditions. Note that we cannot draw firm 
conclusions about the evolution before the transition phase, as this is clearly model dependent. 
If one insists on trusting the lifted exponential potential ([T]) at large 0, then the universe is 
initially in a contracting de Sitter phase, which is of course unstable to kinetic domination. But, 
as mentioned before, there is considerable freedom in specifying the potential in the pre-transition 
phase. In for instance, we present a version of the scenario for which the background solution 
is initially expanding and therefore stable for all times. 



5 Non-Gaussianities and Strong Coupling 

While a phase of rapidly-varying e(t) yields an identical power spectrum as infiation, we will 
see that the degeneracy is broken at the three-point level. Non-gaussianities are strongly scale 
dependent, with the dominant contribution growing as k"^. This implies that our mechanism can 
only generate a finite range of modes within the perturbative regime. A related pathology, also 
discussed below, is that the theory becomes strongly coupled on small scales, invalidating the 
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classical description for these modes. 

These problems all have a common origin: the transition phase with large c is maintained longer 
than needed — a consequence of the simple V{(f)) considered so far. In Sec. [6| we will consider 
altering the pure exponential potential so as to terminate the transition phase before the modes 
with unacceptably large non-linearities are generated, thereby shutting off power on small scales. 
This achieves the desired goal of avoiding strong coupling and large non-gaussianities, while pro- 
viding a range of scale invariant modes on observational scales sufficient to account for microwave 
background and large scale structure measurements. 

5.1 Computing the three-point Amphtude 

For a canonical scalar field with unit sound speed, the exact action to cubic order in ( is given 



where spatial derivatives are contracted with the Euclidean metric 6ij, and x is defined as 



by ini SZl HE] 




+ ^^C'C + f vc ■ Vxv^x + f v^civx)' + 2/(0 

2 2a 4a 




(36) 



(37) 



The last term, proportional to the linearized equations of motion. 





can be absorbed as usual through a field redefinition 



C^C + /(C), 



(39) 



where 



/(C) 



+ 



+ H^^ + + v-^( W(v.cv,C))] 

-^[vc ■ vx - v-^(v^v^-(v.cv,x))] . 



(40) 
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At first order in perturbation theory and in the interaction picture, the three-point function is 

(C(t, k,)at, k,)at, ka)) = -^ f dt'([C(t, ki)C(t, k2)C(t, ka), HUt')]) , (41) 

J — oo 

where i^mt = —L^, up to interactions that are higher-order in the number of fields. As usual we 
expand the quantum field ( in terms of creators and annihilators, 

ak,t)=a{t)a{k) + a{t)aH-k), (42) 
with commutation relations [a(k), a"''(k')] = (2TT)^6^{k — k'). Although the upper limit of inte- 



gration in (41) has been set at t = 0, strictly speaking our approximation a{t) ~ 1 assumed 
throughout breaks down at a time t ~ tcnd/c, as discussed in Sec. [3] We have checked that for 
the modes of interest this makes little difference to the final answer, hence we are justified in 
integrating all the way to t = setting a{t) = 1. To simplify the expressions we set Mpi = 1 for 
the remainder of this section. 



The three-point function receives contributions from each interaction term in (36). The dominant 



contributors, it turns out, are the last two terms in the cubic action, both which are of order 
e'^. The next-to- leading contribution is the rj term. We present explicit calculations of these 
contributions and refer the reader to the Appendix for the rest of the calculation. 

• The contributions: The terms give the combined interaction Hamiltonian 

= -jI d'x (^V^C^C^C + 2CVC^cj . (43) 

Applying the canonical commutation relations, the three-point correlation function (41) in 
this case reduces to 

(C(ki)C(k2)C(k3)).3 = z{27Tf6'{k,+k2 + ks)aMCkMCkM 

^ loo..^' 4 U^^^^^ ^ ^ jnq — d^ + p^^"^- + ' (''^ 

where the small imaginary part at t — t- — oo projects onto the adiabatic vacuum state. 



Substituting the mode functions (18) and using (10) for e(t), we obtain 



(C(kOC(k2)C(k3))e3 = -^^i2Ti)H\k, + k, + k,)K - E ^^^1 + 2^:1^2^:3) 



15 



X 



Im I + ^fc,tend-t.an) ^ dt ' + '^^--^K ^^^ I , (45) 

I • J —oo+ie ' '-end-tranj I 



where K = ki + k2 + k^. 

As usual it is convenient to express the three-point function by factoring out appropriate 
powers of the power spectrum and defining an amphtude A as follows 

(C(ki)C(k2)C(k3)) = (27r)^53(ki + k2 + k3)Pc' , (46) 



where is given by (19). Fortunately, the integrand is a total derivative: 



° 3-2ir(t + tend-tran) ^,i,f = - ^ dt ^ ^ '''''' 



(t + tend-tran)^ J -oo+ie dt \ (t + tend- trail)"' / ^cnd-tran 

(47) 

Substituting ^ for tcnd-tran and focusing on the long wavelength limit -ft'ltcnd-trani ^ 1, 
which is appropriate for the modes of interest, the three-point amplitude is 



\ 



Thus this scales as K"^ /Hi and, as we will see, dominates over all other contributions on 
scales K > \Hq\. 

The f] contribution: The interaction Hamiltonian for this contribution is 

i^int = - jd'x ^e77C'C • (49) 



Substituting (10) and (12) for e(t) and ?7(t), respectively, we obtain the three-point contri- 
bution 



(C(kOC(k2)C(k3)). = {27:f6\\^, + k^ + k3)^|^^Im |n(l + ^k,t 



end— tran ) 



I 1 1 (.^ '^cnd— tran)^ 

X / dt^- 

-oo+ie ~l~ '^end— tran) 



kik]{t + tend-tran) - Kkik2k2,{t + tend-tran)' 

, i i^j 

(50) 
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Performing these various integrals, the corresponding three-point amphtude is given by, in 
the long wavelength (-ft'ltend-trani ^ 1) limit. 



-K K I K 
2 



i ijtj / 



(51) 



This contribution scales as K/\Hq\ and is therefore subdominant relative to (48) on scales 
K>\Ho\. 



The remaining contributions, computed in the Appendix, are all suppressed by 1/c^ relative 



to (48). The full three-point amplitude can be well approximated by (48), at least on scales 
|i^o|: 

^ ^ ^ (^E - E ^^^1 + 2^^^2A:3 j . (52) 

As a check, note that this satisfies Maldacena's "consistency" relation [4S]: in the squeezed limit 
^3 ^ ^1 ~ ^2 = k, we have A ^ 0, consistent with our neglecting departures from scale invariance 
in computing the three-point function. (Although derived in the context of inflation, Maldacena's 
relation applies here because our satisfies its two keys assumptions: single field theory and ( — )• 
constant as A; — )■ |1H]-) Instead our amplitude peaks for equilateral configurations, ki = K/3. 
The shape dependence is qualitatively similar to higher-derivative inflationary models [50j. 

Following standard practice, the three-point amplitude translates into a value for f^'^^' , defined 
at the equilateral configuration |50] : 

/nl ^3 - 144^2- (5dj 

Unlike the power spectrum, the three-point function is thus strongly scale dependent: /nl"'' is 
< 0{1) on the largest scales {K ~ |-f^o|) and grows as K^. This is in stark contrast with the 
small and nearly scale invariant /nl predicted by single-field, slow-roll inflation. The degeneracy 
of our mechanism with inflation established at the two-point level is therefore strongly broken at 
the three-point level. 



The strong scale dependence of (53) implies that perturbation theory breaks down on relatively 
large scales. Specifically, the perturbative expansion parameter is /nlC? with ( ~ 10"^, hence 
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non-linearities dominate for K > 10^/^|iJo|- We will have more to say about this in Sec. 5.3 
In fact, we will see in Sec. 5.4 that on even smaller scales {K > 10^\Hq\) quantum corrections 
dominate the classical result, indicating strong coupling. All of these problems have a common 
origin: the transition phase with large c is maintained longer than needed. A simple solution, 
discussed in Sec. [6} is to alter the pure exponential potential so as to terminate the transition 
phase before these problem emerge. This restores perturbative control in two ways: 1. altering 
the evolution of e(t) suppresses the dominant contribution, thereby expanding the range of 
perturbative modes; 2. terminating the transition phase suppresses ( on smaller scales — the 
spectrum tilts strongly to the red and then flattens out at an exponentially smaller amplitude 
with an acceptable non-gaussianity (/nlC ^ 1) throughout. This leaves us with a finite range 
{\Ho\ ^ K < 10^\Ho\) of scale invariant modes, which is sufficient to account for observations. 



5.2 Horizon Crossing vs Long Wavelength Approximations 



A standard, back-of-the-envelope method for estimating /nl is to compare the cubic and quadratic 
actions for ( at freeze-out [51] : 

^2 



(54) 



/nl ~ 

freeze— out 

This is because non-gaussianities typically peak at horizon-crossing — deep inside the horizon, 
modes are approximately in the vacuum state, whereas far outside the horizon, interactions are 
suppressed by time derivatives and spatial gradients, which are small relative to Hubble in that 
regime. 

The situation in our case is trickier: because the various cubic interactions have coefficients 
such as or e'^ that grow rapidly in time, there is a competition between this growth and the 
derivative suppression. It turns out that for most of the cubic interactions these two effects nearly 
balance out, such that the horizon-crossing approximation provides a good estimate. Consider 
the r] contribution, for concreteness. Since time and spatial derivatives are comparable at horizon 



crossing, we can approximate ( ~ VC ~ k( in evaluating (54) and obtain 



-1 53 
So 



c 



fciti=i 



-1 eiC'C 



fcltl=l 



c 



k\t\=l 



k 



(55) 



where in the last step we have used ?) ~ 2ifQ ^ ~ k^H^ ^ at horizon crossing. This result 



agrees with the parametric dependence of (51). Similarly for all contributions calculated in the 
Appendix. 

The key exceptions are the terms. The vertex increases as ~ during the transition 
phase, and this rapid growth overwhelms the derivative suppression. The three-point contribution 
therefore peaks at late times, well after the modes of interest have frozen out. Indeed, the horizon- 
crossing approximation fails miserably in this case: 



3 53 
^2 



/nl ~ C 



k\t\ = l ^4t74 ' 
fc|i| = l ^ ^0 



(56) 



which greatly underestimates the exact answer ~ k'^/Hl. We can shed further light on this 



contribution by first taking the long wavelength limit of the mode functions (18). Up to an 
irrelevant constant phase, the relevant terms are 



where y = t + tend-tran, and 



a = C,{l + hv + ikV + ...] , (57) 



By inspection, it turns out that the dominant contribution to (44) comes from the imaginary part 
of the integrand: 

+ real part} . (59) 

As expected, this decreases in time due to the derivative suppression. The trouble is that the 
e'^ ~ grows even faster. Indeed, the three-point correlation is 

(C(ki)c(k2)c(k3)).3 = - n . + + ( E - E ^-^1 + ^^^^^^3 

^cnd — trail (^y 



X / — + C.C 

-00+ie y 
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128 n, 



end— trail 



J2kf -^kik'^ + 2kik2h\ . 

(60) 



Substituting the expression ([T]) for tcnd-tran, we obtain the amphtude 



A. 



32H^ 



^k^i -^kik] + 2kik2h 



3 ) 



(61) 



which agrees precisely with (48). In other words, because the y-integral strongly peaks at tcnd- 
the long wavelength approximation reproduces the exact answer. 



tran 5 



5.3 Classical Perturbation Theory 

We have seen that the perturbative expansion parameter /nlC grows larger than unity on small 
scales. This breakdown of (classical) perturbation theory can be seen in other observables, such 
as perturbations in the energy density, Sp/p, in synchronous gauge. (An equivalent discussion 
apphes to the extrinsic curvature perturbation in this gauge, given by ( [35| ), but here we choose to 
focus on Sp/p to parallel the discussion in [S].) With p ^ SHqMpi during the transition phase, 
we have, at linear order, 

5p _ 050 + V^S^ kH 



where we have used (|3| and (31). And since e ~ during the transition phase, this clearly 
peaks at the end of the transition phase: 

rC, (63) 



P 



— fTTid — tran 



in agreement with the parametric dependence in (48). 



The growth in 5p/p at first sight seems to contradict the attractor property established in Sec. |4| 
However, this is an artifact of p being accidentally small: large kinetic and potential energy 
contributions nearly cancel on the background solution, resulting in a comparatively small total 
energy density. If we instead consider second-order contributions, such as 5p'^'^^ ~ 50^, we obtain 



5p(i) 50 + 1/^50 Hq 
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C , (64) 



which clearly becomes increasingly small in time, consistent with the attractor property of our 



solution. Note that evaluating (64) at horizon crossing gives 

5p(2)/5p(i) ^ A;C/|i/o|- Hence the 
perturbation expansion for 5p/p breaks down for k > 10^|ifo|; consistent with the r] contribution 
to the three-point function — see (51). It will be shown below that 10^|ifo| also coincides with 
the onset of strong coupling. 



The authors of |1T] performed a similar analysis in Newtonian gauge and instead found 



6p 
P 



(65) 



Newtonian 



where $ ~ 10~^ is the gravitational potential. If true, then for the values of c of interest this 
would invalidate perturbation theory on all scales. However, this is clearly an artifact of a poor 
gauge choice. As emphasized in [10], Newtonian gauge is ill-suited to study the evolution of 
perturbations in ekpyrotic cosmology, since 50 and $ both diverge as 1/t. In the case of Sp/p, the 
relation between Newtonian and synchronous gauge is (using a{t) ~ Ir^ 



5p 
P 



Newtonian 



^ — ^cnd — tr. 



Sp 

p 

6p 

P 
6p 

P 



sync 



P 
P 

+ 2eC 



H V2 



dt 



H 



sync 



(66) 



sync 



And since $ = at the end of the transition phase [Hj, this is consistent with (65). The large 



contribution in (65 ) is therefore purely a consequence of a breakdown of Newtonian gauge. Similar 



conclusions apply to other quantities in the two gauges, such as 



5.4 Quantum Corrections 

Next we turn to quantum considerations and argue that the pure exponential case studied thus 
far is dominated by quantum effects on small scales. Specifically, we will see that loop corrections 
to the two-point function overwhelm the tree-level contribution, signaling strong coupling. 

^We thank Alex Vikman and Guido D'Amico for discussions on this point. 
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A quick estimate of the magnitude of loop corrections can be obtained by comparing the cubic 
and quadratic action for ( at freeze out, where quantum effects are most important As shown 
in Sec. |5.2[ the dominant contribution at freeze out arises from the r] vertex: 



S3 
S2 



k 



k\t\=l 



(67) 



Thus the theory is strongly coupled for k > 10^|i7o|- 



On yet smaller scales, the stress tensor of quantum fluctuations dominates the background energy 
density, indicating a backreaction problem. This can be estimated by comparing the quadratic 
action 5*2 at horizon crossing with the background action Sq ~ H^Mp^. This gives 



S2 
So 



k\t\=l 



fciii=i 



k' 



(68) 



where in the last step we have used C ~ c\Hq\/M-p\ for the scale invariant modes generated during 
the transition phase. Hence this ratio is also ^ 1 on sufficiently small scales. 



5.5 Summary 

Let us briefly recap the issues uncovered in this Section. We have found that non-gaussianities are 
strongly scale dependent, resulting in a breakdown of classical perturbation theory on small scales. 
The dominant contribution to /nl comes the vertices in the cubic action, which, remarkably, 
peaks at late times, well after the modes of interest have frozen out. As a result, the perturbative 
expansion parameter /nlC becomes larger than unity for 

fc> 10^/2 |^^|_ (gg) 

On smaller scales, loop corrections eventually dominate the two-point function. Specifically, the 
theory is strongly coupled when modes with 

fc>105|iJo| (70) 

are generated. Correspondingly, we have found that the (classical) perturbative expansion for 
5p/p breaks down on those scales. On yet smaller scales, quantum backreaction effects overwhelm 
the background. 
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As mentioned earlier, these problems all result from maintaining the transition phase with large 
c longer than necessary. We will see in the next Section that strong coupling can be avoided 
by altering the potential such that ( strongly tilts to the red on small scales. Suppressing the 
small-scale amplitude in this way in turn allows perturbation theory to be valid on all scales, both 
classically and quantum mechanically. 



6 Weakly- Coupled Model 

The aforementioned small-scale suppression of power can be achieved by generalizing ([T]) to 

V{(j)) = Vo{l- e-'^^*)"^/*^") , (71) 

where c(0) decreases smoothly to 6 ^ c after the transition phase has generated an acceptable 
range of scale invariant fluctuations. It is reasonable to expect that our results will depend on how 
rapidly c(0) decreases and on its asymptotic value b, but should otherwise be insensitive to the 
details of this process. Hence, instead of specifying c(0) it is more convenient to choose a suitable 
e(t) that allows us to proceed analytically. 

We require that e{t) ^ corresponding to c{(j)) = c, from the onset of the transition 

phase until some time tc- Therefore, during the interval tbcg-tran < t < tc, the transition phase 



proceeds as before, and scale invariant modes are generated with amplitude given by (18). This 
standard part of the evolution will be referred to as the scale invariant phase. We assume that e 
is continuous at t = tc and subsequently grows as a power-law: 

where a > 0. (The power-law form is convenient because the C, mode function equation can be 
solved analytically in terms of Hankel functions.) We will refer to this phase as the a phase. 
Meanwhile, since H[t) ~ Hq during the transition phase, we can integrate the relation e(t) = 
-H/H'^ to obtain 

Paradoxically, e(t) increases faster than in the pure exponentially case, which at first sight would 
seem to exacerbate the problems encountered earlier. In fact, this is not so. A faster growth in e 
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can result in a shorter transition phase, which in turn imphes a smaller value of e at the onset of 
the ekpyrotic scaling phase. 



We assume that the a phase ends at a time tg, at which time the universe enters an ekpyrotic 



scaling phase with e = 6^/2. Assuming continuity of e at t^, (72) implies 



2(l+a) 



(74) 



where we have introduced the notation k^ = |tc|~^ and \ks\ = \ts\~^ for future convenience, 
corresponding to the shortest-wavelength modes generated during the scale invariant and a phase, 
respectively. 

During the ekpyrotic phase {t > tg), the Hubble parameter is given by 

m = JTTT^. . for t > , (75) 

" yf^ crunch J 

where tcrunch marks the time of the big crunch. (In reality, we of course envision that the ekpy- 
rotic phase itself terminates before the big crunch and is followed by a bounce to an expanding, 
radiation-dominated phase. In the New Ekpyrotic scenario [21], for instance, a non-singular bounce 
is achieved through a ghost condensate [52J. See [53] for a recent supersymmetric extension of 



this theory.) Matching (73) and (75) at tg gives 



^crunch i o tt ' (76) 

Before turning our attention to perturbations, we note in passing that the modified evolution 
described above can circumvent one of the criticisms raised in [H] , namely that H,H, . . . all 
eventually become super-Planckian in the pure exponential case. Indeed, at the end of the end of 
the ekpyrotic scaling phase, we have Hek-end = c^-f^fk-end/^? which is ^ Mp^ for the values of c 
considered here. Higher derivatives of H are even more singular: 



d"i/ 



dt" 



1 



~ C"+i |i/ek-end| > c'|i7ek-end| > Mpi . (77) 



With the modified evolution, however, the universe eventually matches on to an ekpyrotic scaling 
phase with a much smaller e. All time-derivatives of H will remain sub-Planckian provided that 
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< ^Pi- But b Hiust also be > 1, since an ekpyrotic phase, by definition, corresponds 
to an equation of state parameter larger than unity. In other words, the allowed range is 

i<''^<n^^. (78) 

I -nek— end I 

which can be satisfied for a wide range of parameters. 



6.1 Mode Functions 



Next we solve for the curvature perturbation, tracking its evolution throughout the a phase and 
subsequent ekpyrotic scaling phase. 



a phase: During the intermediate phase in which e evolves as (72), the scale factor remains nearly 



constant, and hence z = a\/2e ~ l/ltl^"*"". The evolution equation (15) therefore reduces to 

^2 (l + a)(2 + a) 



Vk + { k 



Vk = 



for tr < t < t. . 



(79) 



Let us first discuss modes that freeze out during the scale invariant phase (t < tc), ^-e. those 
with k < kc- These modes are already frozen out by the onset of the modified transition phase. 



hence (57) applies just before t = 

CkKkcit < tc) = 



ic\Hn\ ( ^ 1,99 ^ , 

— 1 + -kH^ + -k^t^ 

2V2Mpik^/^ V 2 3 



^0) 



Here we have used \t\ > \tc\ ^ |tend-tran|- By comparing (16) and (79), we note that for a ~ C(l), 
the freeze-out radius | -ff freeze I = \f^T^ only changes by a factor of order unity at tc-, hence modes 



with A; < /cc do not re-enter the freeze-out horizon. Thus, right after t = tc, we can solve (79) in 
the long wavelength limit 



Ck<kcit > tc) — 



1+777^ ^ + 



+ Bu [(-H)3+2« + . . .] 



2(1 + 2a) 

Matching C, and C, aX t = tc allows us to fix and B^. The relevant terms are 



^1) 



k"X 



ic\Ho 



1 + 



PI 



kH^ I fkc 

2(1 + 2a) ^ 3 + 2a v"^ 
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2a 



:-kt) 



3+2a 



(82) 



These modes therefore remain scale invariant throughout the a phase. 



Modes with k > kr, on the other hand, are still in their adiabatic vacuum at the onset of the a 



phase. With this vacuum choice, the mode function solution is 



'Vk>kc 



2Mi 



PI 



a+3/2V ■ 



Using the asymptotic expansion of the Hankel function, the long-wavelength curvature perturba- 
tion = Vk/z on these scales is 



ic\Ho\ 






(1) 



V 



+ 



2(1 + 2a) 22{i+°)(3 + 2a)r2(a + 3/2) 
Hence the spectrum has a strong red tilt for a ~ as desired. 



Ekpyrotic scaling phase: During the ekpyrotic scaling phase {t > tg), the equation of state is 
nearly constant and large, e = 6^/2 ^ 1, and hence the scale factor slowly contracts as power-law. 



a{t) ~ (— t)^/'' . The evolution equation (15) in this case reduces to 

Vk + (k^ - 



Vk = 



for t > tg ■ 



5) 



b (t ^crunch) 

Unlike the scale invariant to a phase transition, the a to ekpyrotic scaling transition typically 



implies a substantial change in the freeze-out horizon. From (85), the freeze-out horizon at the 
onset of the ekpyrotic scaling phase is 

1 



b\t,-t 



crunch 



b\Ho 



(86) 



where we have used (76). On the other hand, for a ~ we can read off from (79) that 

\tg\. (87) 



\t=t7 - 



For our parameter choices, we will see that ^\ /H^ ^L_y-+ = bHots <^ 1, hence some of the 
modes generated during the scale invariant and a phases re-enter the freeze-out at t = t^. We 
must therefore carefully keep track of their evolution. 



The general solution to (85) is 
Vk{t>ts) = Vk{-t) 



^kJ 1 /t,2_8 (^(i^crunch ~ t)) + B^Y^ /b^^s ('^crunch ~ t)) 
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AkJl/2{k{t crunch — t)) + -Bfc^l/2(fc(^cmnch — t)) 



Ak Sin(fc(tcmnch - t)) - Bk COs(A;(tcrunch " t)) 



where the second step follows because 6^1. And since z = av2e 6 is constant in this phase, 
the curvature perturbation is simply given by 



(kit > ts 



^-1 



Ak sin(A;(tcmnch - t)) - B^ cos(A;(tcmnch - t)) 



i9) 



Matching this to (82) and (84), respectively, and using the fact that k < kg for the modes of 



interest, we obtain at late times (A;|tcranch — t\ <^ 1): 



k"\k 



iclHn 



COS 



2k 



ic\Ho\ fk^Y 2^+'^T {a + 3/2) 
AV2Mpi \ k 



cos 



TT 



2k 



for \Hq\ < k < kc ; 

for kc<k <kg. (90) 



Therefore, aside from acquiring an oscillatory factor, ( maintains its original amplitude throughout 
the ekpyrotic scaling phase. The cosine factor results in oscillations in the power spectrum. For 
this effect to be negligible on the largest scales probed by observations, we demand that the cosine 
be approximately constant over the entire scale invariant range. This will be the case if 



b^\Ho\ 



< 1. 



(91) 



On small scales, ( has a strong red tilt, and assumes a minimum amplitude for the shortest- 
wavelength mode {k = kg) generated during the a phase: 



k'/% 



c\Hn 



k I min 



M, 



PI 



(92) 



Finally, on yet even smaller scales, modes with k > kg freeze out during the ekpyrotic scaling 



phase, and as usual have a strong blue tilt. Imposing the adiabatic vacuum choice in (88) fixes 
Ak,Bkr^ l/V2kMpi, hence 



k 



for k > k. . 



(93) 



bMpi 

This growth is cut off once the ekpyrotic phase terminates, which occurs well before the amplitude 
reaches unity. 
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6.2 Avoiding Strong Coupling 



The strong red tilt on intermediate scales generated during the a phase can cure the strong coupling 



problem encountered in Sec. 5.4 The dominant contribution to S'3/S'2 at horizon crossing is given 

k 



as before by (67) 



^2 



fc|i|=i i-'-'ol 

On the largest scales [k < kc), ( is scale invariant and 
{k > kc), however, we have 

^3 
^2 



(94) 



10 ^ as before. On smaller scales 



k 



k\t\ = l 



< 



(95) 



For a > 1, in particular, the theory becomes increasingly weakly coupled on small scales. Hence, 
provided that the range of scale invariant modes satisfies 



kc < 10' li^ol 



(96) 



then for a > 1 quantum corrections are under control on all scales. This is our main constraint 
on the allowed range of scale invariant modes. 

The a phase also allows us to circumvent the quantum backreaction problem of the pure expo- 



nential case. As in (68), the backreaction is largest on small scales, hence we focus on the modes 



generated during the a phase: 



52 

5-0 



A:|t| = l 



HI 



U2U2 / u 

hj rv^ I r\j 



k\t\=l 



2{l+a) 



(97) 



where in the last step we have substituted (90). Although the parametric dependence is identical 



to (68), the upshot of the a phase is that it limits the range of modes generated. Backreaction 



peaks at k = kg and is under control provided that 

kl<\Ho\Mpi. 



(98) 



6.3 Non-Gaussianities 

The modified evolution for e(t) should have a dramatic impact on the three-point function. Indeed, 
recall that the dominant e'^ contribution peaked at late times, which is precisely what has been 
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altered with the a phase. 



To see how non-hnearities can be tamed, let us focus on the scale invariant modes (|-ffo| < ^ < ^c)- 
During the a phase, their evolution is described by (82). As in Sec. 5.1 , the dominant contribution 



to (44) comes from the imaginary part of the integrand. Substituting (|82|), we obtain 



k\ 



-%)c\Eq\ 



/ 2y2MpiA; 



3/2 



2i 



X 



l + 2a' 



-K^ y^k\-^ kik] + 2kik2ki |tcr^"(-t)^+^" + real part \ . (99) 

While this is suppressed for small t, the integral is once again overwhelmed by the growth in the 
vertex: ~ l/t^*^^"*""). The three-point function is 



{C(k,)C(k2)C(k3)), 



Aa 



X 



1 + 2a 



dt 



t3- 



Aa 



+ C.C. 



(27r)^5nE.k. 



2j.2o 



,2(1+2q) 



128n,A^| (l + 2a)2t 
Rewriting this in terms of kf. = l/|tc| and kg = ^/\ts\, we find the amplitude 



A, 



9M^, klK^ 



with corresponding /^l"'' parameter: 



^ - ^ + 2kik2k^ 



, X 2(l+2a) 

kf 



(101) 
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'NL 



(102) 



Note that, remarkably, if the a phase is maintained long enough to the point where the approx- 
imation H ^ Hq breaks down, that is, if tg is chosen such that \Hq\ ~ 2c~^/cc(l + 2a)(tcAs)^'''^", 
then (102) exactly matches our earlier result for the pure exponential case: /^l ~ k'^/^o- There- 
fore, by terminating the a phase at an earlier time, we can suppress f^l^, as desired. Specifically, 
we demand that /nlC < 1- For the scale invariant modes, \Hq\ < k < kc, the amplitude is 
( ~ 10~^, and the condition is most stringent at A; = kc'- 



(/nlC) 



f-^ \k=kc 



■^vi \k 
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2(l+2o) 



where in the last step we have used (21) 



It turns out that a similar calculation for the kg < k < kc modes, obtained by substituting the 
mode function (84) into the three-p 

/c^/^C ~ in this case — see (90) — , it follows that (/nlC) 



mode function (84) into the three-point amplitude, yields an identical expression for /^^ 



^3 ~ fc^ " peaks at k 



Since 
kr a a > 2, 



and hence is automatically < 1 when (103) is satisfied. For simplicity, we will therefore impose 



a>2. 



(104) 



The remaining contributions to the three-point function can be estimated by the horizon-crossing 
approximation, as in Sec. 5.2, Since r/ ~ 1/t during both the scale invariant and a phases, we have 
/nl ~ k/lHol on all scales. As in the strong coupling discussion of Sec. 6.2, /nlC peaks at k = kc 



for a > 1, and is < 1 provided (96) is satisfied. As before, the contributions are subdominant 
and can be neglected. 



6.4 Summary of Constraints 

To be phenomenologically viable, the generalized model described above must satisfy the following 
list of constraints: 

1. Correct large scale amplitude: The amplitude of the long- wavelength modes {k < kc) 



should match observations of the large-scale power spectrum. From (21) we have 



lo-^ (105) 



PI 



2. Scale invariant modes must match the observable range: The modes generated during 
the scale invariant phase are on scales smaller than |ifo|^^; hence the comoving scale l-f^ol^^ 
must encompass the entire observable universe. This leads to the upper bound on \Hq\ given 



by (23), where recall that -ffek-end, the Hubble parameter at the end of the ekpyrotic scaling phase, 
is assumed comparable in magnitude to the Hubble parameter at the onset of the expanding, 
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radiation-dominated phase: |ifek- 



end I 



^reheat /^Pi- ^oi simplicity, we will assume that the 



bound (23) is saturated, thereby fixing \Hq\ in terms of the reheating scale: 



reheat 



(106) 



3. Avoiding strong coupling: As discussed in Sec. 6.2, loop corrections are small if the range 
of scale invariant modes is restricted to kc < 10^ \Ho\ — see (96). (Another necessary condition 



is a > 1, but this now follows from (104).) Since a factor of 10^ is just enough to be consistent 
with microwave background and large scale structure observations, we fix k^ to nearly saturate 
this bound 

kc ~ 10^ \Ho\ . (107) 



Meanwhile, quantum backreaction is under control provided that kg < |ifo|^Pi 



Using (106) and (107), we can rewrite this as 



k. 



> 10" 



Tr 



reheat 



M, 



PI 



see (98). 



(10^ 



4. Scale invariance and sub-Planckian curvatures: The equation of state during the ekpy- 
rotic phase, e = b'^/2, is constrained by two requirements. To avoid reaching super-Planckian 



curvatures by the end of the ekpyrotic phase, 6^ is bounded from above through (78). Given 
our assumption that |/7ck-end| sets the reheating scale, the upper bound implies 6^ < ^pi/T^^eheaf 



Substituting (74) and using the relations (105) and (106), we can rewrite this inequality as 



> 



10 



-20 



Tr 



reheat 



M, 



PI 



1 

l + Q 



(109) 



The equation of state is also bounded from below by demanding that the scale invariant spectrum 



is not appreciably distorted by the ekpyrotic scaling phase. From (91) and (107), this requires 
6^ > 10^. Once again substituting (74) etc., we obtain 



< 



10 



-45 



-^reheat 
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2(l-|-a) 



:iio) 



5. Validity of classical perturbation theory: We have seen in Sec. |6.3|that, provided a > 2 



the tightest constraint from non-gaussianities comes from the contribution at k 



kc- Us- 



ing (106) and (107), the inequality (103) can be rewritten as a constraint on kc/ks'. 

kr- 



k. 



> 



10 



-85 



T \ ' 
-L reheat i 

Mpi J 



2{l + 2a) 



:iii) 



This is generally a more stringent lower bound than either (108) or (109) 



6.5 Working Examples 



To summarize, once we choose a reheating scale Treheat, the scale of the potential \Hq\ (or equiv- 



alently, Vq) is fixed by (106). In turn, the exponent c characterizing the scale invariant phase is 



fixed by the large scale normalization (105), while the comoving scale kc marking the onset of the 



a phase is determined by (107). Only two parameters remain to be specified: a and kc/kg. For a 



given a > 2, we will check that kc/kg satisfies the inequalities (108) — (111). 



High-Scale Example: Consider GUT-scale reheating, T^eheat = 10"^ Mpi = 10^^ GeV. From ( |105| ) 
and (106) this fixes Vq^^ = 10 GeV and c = 10^^. Choosing a = 5, the upper bound (110) re- 



duces to kc/kg < 5 X 10 ^. Meanwhile, among the lower bounds, (111) is the most stringent: 
kc/kg ^ 4 X 10~^. Dropping factors of order unity, we therefore impose 

Kg 

Hence ~ 12 e-folds of modes are generated during the a phase in this case. 



(112) 



Intermediate-Scale Example: Consider reheating at an intermediate scale, Treheat = 10 ^ Mpi 



10^ GeV, corresponding to Vq^^ = 10 ^ GeV and c = 10^*^. Choosing a = 3, the upper bound (110) 



reduces to kc/kg < 10 ^, while the tightest lower bound is again given by (111): kc/kg > 8 x 10 ^. 
Hence in this case we impose 

^ lo-^ (113) 



kc 



corresponding to 18 e-folds of a-phase modes. 
Low-Scale Example: Consider electroweak-scale reheating, Treheat 



10 



-15 



corresponding to Vq^^ = 10 ^ GeV and c 



10^0. With a 
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Mpi = 10^ GeV, 
3, we find that the inequalities on 



kc/kg are satisfied for 

^ = 10-1°, (114) 
kg 

whicfi amounts to ~ 23 e-folds of a-phase modes. 



7 Conclusions 

In tliis paper, we have explored tlie adiabatic ekpyrotic meclianism proposed recently to generate 
a scale invariant spectrum within an attractor background. At the level of the power spectrum, 
the adiabatic mechanism is dual to inflation — the equation governing C, and its growing mode 
solution are identical to inflationary cosmology. 

As we have seen, however, the duality is broken by the three-point correlation function. Unlike the 
nearly gaussian spectrum of inflation, the rapidly- varying equation of state e ~ 1 /r^ characteristic 
of the adiabatic ekpyrotic phase results in large non-gaussianities on small scales. For the simplest 
exponential potential ([T]), the most dangerous contribution comes from 0{e^) terms in the cubic 
action, which are subdominant in the inflationary case. The rapid growth in these vertices gives 
three-point contributions that peak at late times when e ^ 1. At the same time, loop corrections 
dominate the tree-level computation on small scales. 

This strong coupling and perturbative breakdown both trace back to the fact that the transition 
phase with large c is maintained longer than necessary. As we have seen, these pathologies can 
be avoided by considering more general potentials where the exponent decreases smoothly from c 
to a much smaller value 6 <^ c once a suitable range of scale invariant modes has been generated. 
This suppresses power on small scales, and thereby restores the validity of perturbation theory on 
all scales. We have shown that the resulting range of nearly scale invariant and gaussian modes 
can span at most a factor of 10^ in k space, which is enough to account for microwave background 
and large scale structure observations. 

We are currently generalizing the scenario to the case of time-dependent sound speed Cs(r), as 
expected in non-canonical scalar field theories, with the hope that this can alleviate the issue 
of non-linearities. As shown in [31], there is much more freedom in generating scale invariant 
perturbations in this case: for any background a(r) there exists in principle a suitable Cs(r) such 
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that ( acquires a scale invariant spectrum. It will be interesting to see if Cs{t) can also tame the 
growth in the three-point function. 
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8 Appendix: Calculation of three-point Function 

In this Appendix we complete the calculation of Sec. [5] and compute all remaining contributions to 
the three-point function for the lifted exponential potential ([T]). As in Sec. [sj we assume Mpi = 1 



throughout. The contributions listed below refer to the cubic action (36) 



The (('^ contribution: This interaction Hamiltonian in this case is i^int = — J d^xe^CC^- 



Following similar steps as in Sec. 5A, we obtain 

.2 7t2\ 3 



X 



(C(ki)C(k2)C(k3))^^. = ^(27r)353(ki+k2 + k3)^(^^^ n(l+^Mcnd-tran) 

3 3 j 

J2 e^t) (t + tend-tran)' [k^k^c'"'' - ik^k^k^it + W-tran)fcifc,e*^*] + C.C. 

J -oo+ie 

^ -(27r)353(ki + k2 + k3) ^/^\3 V k^k] log(X|tend-tran|) , (US) 

where in the last step we have assumed -ft'|tend-tran| ^ 1, appropriate for the modes of 
interest, and where we have used the identity 

K 5^ k}k] - k.k^k^ Yl ^i^j = Yl ^i^^ ■ (^^6) 

i<j i<j iy^j 

The corresponding amplitude is 
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i^J 



This contribution is suppressed by 1/c^ relative to the dominant, amphtude, given by (48). 



The C(VC)^ contribution: The interaction Hamiltonian in this case is i^mt = — J d^x e^C(VC)^, 
with corresponding three-point function 



(C(kOC(k2)C(k3))c(vc)^ = -i^^r^H^i + k2 + ks)^^^ Im |n(l + ^k,t 

— — 1 - iKt - 2J kikj{t + tcnd-tran)^ + ikik2k3{t + tend 



end— trail J 



^3 
tran I 



(118) 

where we have used the identity 

ki-k2+k2-k3 + ki-k3 = -^5^A;2. (119) 

i 

Performing the integrals in the -f^ltend-tranl ^ 1 limit, we obtain the amplitude 

-^C(VC)^ ~ -J^^ Yl E ^1 log(^l^end-tranl) • (120) 
i j 

This is also suppressed by 1/c^ relative to the contribution. 



The CVCVx contribution: From the definition of x in (37), we find Hi^t = 2e^ / d^x CVC^C- 
The three-point contribution is therefore given by 



(C(ki)C(k2)C(k3))cVfVx = (2^)'^'(kl + k2 + Im |n(l + zA:,tend-tran) 

X / -— (k2 ■ k3A;^[l - zk2{t + t,nd-tranW' + perms.) \ . (121) 

J — oo+ie ~r ''cnd-tranj J 

Performing the integrals and using the identities 



+ perms. = + 2A;,fc2A^3 + E ^'^1 ' (^22) 
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we find the following leading piece for -ft'] tend- tran I ^ 1 

%CVX = ^ ( E ^'-\Y. + ]^ E ^'^1 + 4^1^2fc3 ) l0g(ir| W-tranI) • (123) 



Again there is a l/c^ suppression compared to the dominant contribution. 

Field Redefinition: The three-point function also receives contributions from the field 



definition (39): 

^2^^^^ ■ ~ v-2(v^v^-(v.cv,x))] . (124) 

Most of these terms involve time and/or spatial derivatives, and hence give negligible con- 
tribution deep in the ekpyrotic scaling phase, when the modes of interest are well outside 
the Hubble radius. The one possible exception is the rjC,"^ term, but this contribution is also 
suppressed deep in the ekpyrotic scaling phase, since e ~ 0^/2 is approximately constant and 
hence r/ — )■ 0. We can therefore safely ignore the contributions from the field redefinition. 



8.1 Summary 



The non-vanishing contributions to the three-point amplitude, given in (48), (51), (117), (120) 



and (123), are 
A. 



V 



A-- - 



n K K 

8 I y 



2cm, 



trail 



E E log(i^|te„d^tran|) 



^ E - I E + E ^^'^1 + 4^1^2fc3 ) log(irltend-tranl) • (125) 



2cmi 
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